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Effects  of  an  Oscillating  Torsional  Couple 
on  the  Contact  Surfaces  of  Elastic  Spheres 


Introduction 

The  problem  of  two  elastic  bodies  in  Hertz  contact  subjected  to  a 
monotonically  increasing  torsional  couple  applied  about  the  axis  of  sym- 
metry of  the  system  was  first  considered  by  Mindlin^  who  found  that, 
in  the  absence  of  slip  on  the  interface,  the  tangential  traction  rises  to 

infinity  on  the  elliptic  boundary  of  contact.  The  torsional  contact  of 

(2) 

two  like  spheres  was  reexamined  by  Lubkin  who  made  allowance  for  the 
occurrence  of  slip. 

The  present  paper  extends  the  Lubkin  solution  to  the  case  in  which 
the  normal  force  is  held  constant  while  the  torsional  couple,  after  having 
reached  a certain  value,  decreases,  and  to  the  case  in  which  the  torsional 
couple  oscillates  between  fixed  amplitudes.  In  the  latter,  the  curve  re- 
lating twisting  moment  to  angle  of  twist  forms  a closed  loop,  the  area  of 
which  represents  the  frictional  energy  dissipation  per  cycle.  For  small 
amplitudes,  this  energy  loss  is  found  to  vary  as  the  cube  of  the  twisting 
moment.  The  solution' is,  strictly  speaking,  valid  only  for  small  amplitudes 
of  the  applied  moment,  because  of  an  approximation  made  at  the  start.  Prac- 
tically, however,  it  is  good  over  a large  portion  of  the  permissible  range 
of  moment-normal  force  ratios. 


(1)  R.  D.  Mindlin,  "Compliance  of  Elastic  Bodies  in  Contact,"  J.  Appl. 
Mech.,  Vol.  16,  1949,  pp.  259-268. 

(2)  J.  L.  Lubkin,  "The  Torsion  of  Elastic  Spheres  in  Contact,"  J.  Appl. 
Mech.,  Vol.  18,  1951,  pp.  183-187. 
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An  approximate  expression  is  obtained  for  the  torsional  compliance 
which  depends  upon  the  past  history  of  loading.  Ihe  initial  value  of  the 
compliance  during  the  first  loading  process  is  found  to  be  the  same  as 
given  in  reference  (1). 

Summary  of  Previous  Work 

(3) 

The  Hertz  problem  for  two  like  spheres  predicts  a plane , circular 
contact  surface  of  radius 

a = (QNR)  / (i) 

where  N is  the  normal  force,  R the  radius  of  the  spheres,  and  Q - 3 (f~  , 

ir.  which  V and  £ are  Poisson’s  ratio  and  Young’s  modulus,  respectively,  cf 

the  material. 

/ 

The  distribution  of  normal  traction  on  the  contact  surface  is  given  by 

cr  = (2) 

zwor  f 

where  yo  is  the  radial  distance  from  the  center  of  the  contact  surface. 

If,  now,  an  additional  moment  ( M ) is  applied  about  the  axis  of 

normal  contact,  and  if  it  is  assumed  that  no  slip  occurs,  symmetry  considera- 
tions lead  to  the  conclusion  that  no  normal  component  of  traction  is  induced 
and  that  the  entire  contact  surface  rotates  as  a rigid  body  with  respect  to 
a distant  point  of  one  of  the  spheres.  The  situation  is  thus  reduced  to  a 
mixed  boundary  value  problem  in  elasticity.  With  a system  of  cylindrical 
coordinates  whose  origin  is  at  the  center  of  the  circle  of  contact  and  the 

(3)  S.  Timoshenko  and  J.  N.  Goodier,  Theory  of  Elasticity,  (McGraw-Hill 
Book  Company,  New  York,  1951),  p.  372. 
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Z-axis  along  the  axis  of  symmetry  of  the  two  spheres,  the  boundary  condi- 
tions are:  on  the  contact  surface,  the  normal  component  of  traction 

(zero)  and  the  tangential  component  of  displacement  (proportional  to 
the  distance  from  the  origin)  are  given.  On  the  remainder  of  the  surface 
of  the  sphere,  approximated  as  plane,  the  three  components  of  traction 

* (all  zero)  are  given.  Solution  of  this  problem ^ gives 
the  shearing  stress  Tz0  on  the  contact  surface, 


r m 


L 


4v'a* 


(3) 


and  the  torsional  compliance  of  a single  body 


(4) 


where  ft  is  the  angle  of  rotation  of  the  contact  surface  with  respect  to  a 
distant  point  in  the  body  H ^ O j and  is  the  shear  modulus  of  the 
material. 

/ 

Since  t'  becomes  infinite  on  the  boundary  of  the  contact  circle,  slip 
must  result  from  an  applied  twisting  moment,  no  matter  how  small.  It  may  be 
expected  to  start  at  the  edge  of  the  contact  surface  where  the  singularity 
occurs,  and,  because  of  the  radial  symmetry  of  T , to  progress  radially 
inward  on  an  annulus.  Further,  it  is  assumed  that,  on  the  annulus  of  slip, 

T 9 y <r  , where  f is  a constant  coefficient  of  friction,  and  that  it 

does  not  exceed  this  value  elsewhere.  On  account  of  symmetry,  the  rigid 
rotation  remains  unaffected  for  the  region  over  which  no  slip  has  occurred 
(the  "adhered”  portion).  This  gives  rise  to  another  mixed  boundary  value 
problem  in  elasticity.  The  normal  component  of  traction  (zero)  and  the 
tangential  component  of  displacement  (proportional  to  the  distance  from 


the  origin)  are  given  on  the  adhered  portion  and  the  traction  is  given  ever 
the  renainder  of  the  boundary  ( TMr  -<?£--  C on  the  annulus 

of  slip  and  all  components  of  traction  zero  outside).  Solution  of  this 

(2 ) 

problea  yields  the  shear  stress 


r’  ^ (aKr'y/'< 

v v>  5 

ZnoF  1 <■  rrL  ..  .7) 

- :/K(*)£(*rv)/f  5 

where  C is  the  inner  radius  of  the  annulus,  <'s  c/a  , K - ( / - 
Ft*',  *>)  and  £(k'ttp)  are  incomplete  elllptio  integrals  of  the  first  and 
seecnd'klnd,  re  s pe  e t ive ly , - of  - modulus  xj  and  amplitude 

r r 

and  D(k)  1b  a tabulated^  ooopleta  elliptic  integral,  of  rcdulua  <, 

ID  * ( K - £ V*.1,  (6) 

Here,  /f^  and  tE  are,  respectively,  the  complete  elliptic  integrals  of 
the  first  and  second  kind,  of  modulus  K.  . A eroes-section  of  the  distribu- 
tion of  V is  illustrated  by  curve  a A 0 in  Figure  1. 

The  relation  between  the  applied  moment  and  the  inner  radius  of  slip 
is  given  by  the  condition  of  equilibrium 


a 

M ~ 27T  /Z7V/4 5 


(4)  E.  Jahnke  and  En.de , Tables  of  Functions,  (Dover  Publications, 

Hew  York,  1945),  pp.  V and  83. 
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Henoe . 


f *_«■  / 3 
4 7T  7 ^ 


/ ^ + k'k'ClIK  + (4  «'v-  *)c]  - 3*.(K 

rfi 

4 3(C*"/  K / 

L 0 (/-*t\ 4^V)Vi 

*A  , 

• ID  f ' (*"' ***•»<)  dot  I I 

0 (,-*'V~SO'A  J > 


Further,  the  relation  between  the  angle  of  twist  and  the  inner  radius  of 
slip  is 

A * /c*  fD  (9 ) 

fyu*1 

Equations  (8)  and  (9)  express,  implicitly,  the  relation  between  the  applied 
twisting  moment  and  the  resulting  angle  of  twist,  'This  is  represented 
graphically  by  curve  OPF  in  Figure  2,  The  value  of  the  moment  at  which 
slip  ooc\irs  over  the  entire  contact  surface,  i.e,,  sliding  is  initiated,  is 
ftjjNa,  * 3rr//6  . This  corresponds  to  an  infinite  value  of  the  angle 


of  twist. 


///Va  « / 


In  what  follows  it  will  be  necessary  to  have  explicit  and  reversible 
relations  between  each  pair  of  the  three  quantities:  inner  radius  of  the 

slip  annulus,  twisting  moment  and  angle  of  twist.  This  is  manifestly 
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impossible  for  the  exact  solution  [Eqs.  (2)  and  (9)j«  us>  therefore, 

confine  cur  efforts  to  the  range  of  small  M ’ i.e»»  : c/a  * / • 

Because  of  the  divergence  of  the  two  integrals  in  Eq.  (8)  for  K>  = / 

all  efforts  to  express  M of  Eq.  (8)  as  a polynomial  in  < have  been  un- 
successful, For  this  reason  the  equation  for  M has  been  obtained  in  a 

different  form.  3tarting  with  the  following  expression  for  V > 

f e (f.i. 

s r > r 


T ■a 


(10) 

/(<zx'f>1) 
a 1 K^a***** /{a1-/*  V]  ( I ~ 


dt<j 


p &c 

inserting  it  into. the  condition  of  equilibrium,  Eq.  (?)»  reversing  the  orSer 
~of~ integration,  and  noting  that ^ ^ 

*k  _ 

f -£2JL—  ^ * /D 

« l/  / - ACV^V 


r/t  - 


<+■ 


K 


we  find 

j2L 

//v* 


•b 


k,  k'*c~  £ 

n 

rk 


(n) 


(5)  In  reference  (2),  insert  Eq.  (47)  into  Eq.  (36).  The  a'  is  here 
replaced  by  c . 

(6)  L.  Fotin,  Formulas  et  Tables  K'Wriques,  (Gauthier-Villars , Paris, 
1925),  pp.  712,713. 


This  expression  is  easily  expanded  in  a power  series  in  K.  with  the 


result 


a ~ „*/,•  _ £*.*-  -—k*  - . , ■) 

JnZ  ~ * < ' * i*  > 


Table  I.  compares  the  values  of  M/fNc,  , calculated  b7  retaining  two  and 
three  terms  in  the  expansion,  2a.  '12 ))  'with  the  exact  values  tabulated  in 
reference  (2). 

Table  I. 


It^e/cv  KLei-kll 


M /fN<* 


Exact 

3 -term 
Approx, 

% Error 

2-term 

Approx, 

0 

0 . 

0 

0 

0,0938 

0.09392 

0.1 

0.09393 

0.1764 

0,17636 

0 

0.17646 

0.2482 

0.24829 

0 

0.24862 

0.3105 

0.31067 

0.1 

0,31140 

0.4082 

0.  U039 

0,5 

0,41246 

0.4805 

0.48230 

0.4 

0.48640 

0.5286 

0.53247 

0.7 

0.53906 

0.5592 

0.56613 

1.2 

0.57540 

0.5769 

0.58768 

1.9 

0,59946 

0.5855 

0,60058 

2.6 

0,61440 

0.5890 

0.60937 

3.4 

0.62  5 00 

% Error 


WOHO 


A similar  expansion  of  in  Eq.  (9)  yields 

bifH  ar  1 £(t  * £ *.*•  + ■ • ) 

jH  /6  V o 6f  / 


(13 


but  its  use  must  be  ocnfined  to  small  K <,  since  the  expression  in  Eq.  (9) 
diverges  for  /C  : / . Exaot  and  approximate  values  of  are 

compared  in  Table  II. 


Table  II. 


/ - K,v 

- Esaat"' 

.3  •term 
Approx. 

t Error 

2-term 

Approx, 

■ ?-Error 

1' 

0 

O '." ' 

" O' 

0 

0 

0 

0.95 

0.0975 

0.0190 

0.01899 

0 

0,01895 

0.3 

0.9 

0.19 

0.0385 

0.03847 

0.1 

0.03817 

1.0 

0.85 

0.2775 

0.0586 

0*05838 

0.4 

0,05744 

2.0 

0.8 

0,36 

0.0794 

0.07866 

0.9 

0,07661 

3.5 

0.7 

0.51 

0.1234 

0.11975 

2.9 

0.11392 

mBSm 

0.6 

0,64 

0.1716 

0.16527 

3.7 

0,15395 

Bl 

In  view  of  the  great  complications  arising  from  employing  three  or 
more  terms  in  the  expansions,  Eqs,  (12)  and  (13),  ve  confine  o'urselves  to 
two- term  approximations.  A short  computation  yields  the  following  relations: 


j 

! 


Equation  (17 ),  which  represents  the  explicit  moment- twist  relation, 
lead  directly  to  the  torsional  compliance 

3 r.  ( . *._M_  r . 1 

r ' /i/ja‘L  2 ( / ' * / ~ J (1S) 

It -is  seen  that  it9  initial  value  is  the  same  as  in  Eq*  (4), 


|k 

Suppose,  now,  that,  after  having  reached  a value  , suoh  that 
0 4,  < 2sTr /tt>  , the  twisting  moment  ft)  is  reduced.  If  slip 

were  prevented,  the  tangential  component  of  traction  V would  tend  to 
negative  infinity  on  the  boundary  « a • This  conclusion  is  reaohed  from 
the  solution  of  the  appropriate  boundary  value  problem  in  elasticity,  as 
described  in  the  asctior.  entitled  ''Summary  of  Previous  Work".  Hence,  slip, 
opposite  in  sense  to  the  initial  slip,  is  presumed  to  start  at  ^ •=  a,  and 
penetrate  to  a radius  b , assumed,  temporarily,  to  be  of  magnitude  c^k^CL. 
On  the  annulus  ^ yia  } "T"  r - j cr  , so  that  the  change  of  tangential 

(7)  The  procedure  in  this  and  the  following  sections  was  first  employed 
in  connection  with  tangential  contact  of  two  spheres.  See 
R,  0.  Mindlin,  W.  F.  Mason,  T,  F.  Osmer  and  H.  Ceresiewicz , "Effects  of  an 
Oscillating  Tangential  Force  on  the  Contact  Surfaces  of  Elastic  Spheres," 

Proc.  1st  U,  S.  Nat’l.  Cong.  Appl.  Yach.  (1951),  pp.  203-208. 
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traction  is  - Ljtr . Since  no  additional  slip  occurs  on  the  surface 
the  chance  of  displacement  in  that  region  must  be  that  of  a rigid  body  rota- 
tion, Thus,  the  change  due  to  reduction  of  A7  presents  a boundary  value 
problem  in  elasticity  identical  in  form  with  that  solved  in  reference  (2), 
Hence,  by  analogy  with  So.  ( 5 ) » the  change  in  traction  is 

TC  * “ 

r 2i£L  rf1*  f i * DMF  <«/,»,)  (19> 

7 rev*  i i 7r  l 

9,)]  j , 


where  , ftj  and  are  obtained  frees  K , and  ^ , respectively,  by 
replacing  C by  k . This  distribution  is  illustrated  by  ourve  In 

Figure  1. 

The  resultant  traction  accompanying  a reduction  in  M is  obtained  by 
adding  the  initial  traction,  Eq«  (5),  and  the  charge,  Sq.  (19),  with  the 
result  t 


jy  (O  ± <X. 


t„  r - iddj  ('  + $At)  , Cif>  «=*  (20) 

2 TT  A 

Va  - - + £(t\  -*<)],  r-c 


where  A ^ stands  for  the  bracketed  expression  in  Sq.  (19)  and  Ac  for  a 
similar  expression  in  Eq.  (5).  The  traction  profile  is  illustrated  by  curve 
&DEO  in  Figure  1, 

The  condition  of  equilibrium,  Eq.  (7),  yields,  upon  manipulation  des- 
cribed in  the  section  entitled  ''Approximate  Formulas'1, 


(21) 


M - 


no*  [, . ( 


i + J i>  y<aL 

* 


Noting,  from  Eq.  (14),  that  the  first  term  on  the  right-hand  side  of  Eq.  (21) 
is  the  -initial  torsional  moment  M*  , the  inner  radius  of  the  annulus  of 
counter-slip  is 


3 /j*-  £2  ) 
Z 2/Na  / 


(2 


2) 


Thus,  as  long  as  - M*  6 M ^ , the  assumption  c 6 is  valid. 

a* 

When  At  s - M , i,e,,  when  the  twisting  moment  is  fully  reversed,  b a c, 
i.e,,  counter-slip  has  penetrated  to  the  depth  of  initial  slip.  At  the  same 
time,  Eos.  (20)  reduoe  to  Eqs.  (5)  with  signs  reversed.  The  traction  is  then 
distributed  Just  as  the  initial  traction  at  A1  * A1*waa,  but  'with  opposite 
sense,  The  sit'uation  is  the  same  as  if  no  positive  Ad  had  ever  been  applied, 
'out  only  a negative  M of  magnitude  Ad4, 

The  associated  rotation  of  the  adhered  portion  is  found  by  a similar 
process  of  superposition.  First,  the  change  in  angle  of  twist  is  obtained  by 

multiplying  Eq.  (16)  by  -2  and  replacing  c by  b as  given  in  Eq,  (22).  The 

initial  twist  is  given  by  Eq.  (16)  in  which  the  value  of  c is  found  from 
Eq,  (15)  wherein  A1  is  replaced  by  Ad*  . The  resultant  angle  of  twist, 

given  by  the  sum  of  the  two,  is  given  by 


- 3 

2 ZjNcx 


r- 


['*(• 


(23) 

3 _Ad, 
y/Va 
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The  relation  between  twisting-moment  and  angle 


of  twiit  is  illustrated  'ey  the 


full  lines  OP RS  in  Figure  2. 

The  compliance  for  unloading  is 


/m"-m  yl/l 

ZfN<k  J 


] 


(24) 


Thus,  the  initial  compliance  on  unloading  [ M - FA*  in  Eq,  (24)]  is  the 
sane  as  the  initial  compliance  on  loading  ] M -0  in  14..  (18)],  That  is, 
in  Figure  2,  the  slope  of  PR  at  P is  the  same  as  the  slope  of  OP  at  O. 

When  PA  has  been  reduced  froc  PA*  to  zero,  there  is  a permanent  set 
given  by  OR  in  Figure  4>  the  magnitude  of  which  is  obtained  by  setting  PA*0 
in  Eq.  (23).  The  accompanying  traction  is  not  zero,  but  is  a self-equilibrating 
distribution  obtained  by  setting 


Jl 

* fNCK 


I 


in  Eqs.  (20). 

when  M has  been  reduced  to  “P'1*,  the  twist  has  reached  the  negative 
of  the  twist  at  MsM*(i,e.,  the  abscissa  of  5 in  Figure  2 is  the  negative 
of  the  abscissa  of  P ) and  the  compliance  is  identical  with  that  of  curve  OP 
at  PA  s FA*  < Hence,  the  unloading  curve  PR 5 is  tangent,  at  5 i to  the 
central  perversion,  OS  , of  the  loading  curve  OP > 


Oscillating,  Torsional  M spent 

Ln  the  preceding  section,  the  entire  situation  at  (A  = -A1*is  identical 

0 

with  that  at  fA  - PA  , except  for  reversal  of  sign,  Hence,  a subsequent  in- 
crease of  fA  from  ~ PA*  to  FA*  will  be  accompanied  by  the  same  events  as 
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occurred  during  the  decrease  from  M to  - M , except  for  reversal  of 
sign.  Thus,  in  starting  along  SU  , Figure  3,  the  compliance  at  5 is  the 
same  as  the  compliance  of  PR  at  P . Slip  again  starts  at  = « in  the 
sane  sense  as  occurred  along  path  OP  . At  an  intermediate  point  of  SU  the 
traction  will  be  like  &DEO  , Figure  1,  with  sign  reversed.  When  the  twist- 
ing moment  once  more  reaches  M*  , k will  again  have  penetrated  to  C and  the 
traction  will  oe  exactly  a AO  , Figure  1. 

The  twist  along  path  SUP  is 


“A  ^ M ) 


Henoe,  when  M s PI*  , the  terminal  points  P of  SUP  and  P of  OP  (Figure  3 
are  identical,  aS  nay  be  seen  by  comparing  So.  (25)  with  So.  (17)  for  that  value 
of  M . 

Subsequent  diminution  of  /H  will  then  produce  a repetition  of  the 
events  accompanying  the  first  diminution. 

It  is  now  evident  that,  during  oscillation  of  PI  between  , with 

N maintained  constant,  the  torque-twist  curve  traverses  the  loop  P R S U P , 
Figure  3. 

The  area  enclosed  in  the  loop  gives  the  frictional  energy  loss  per 

cycle : 

F = L (k  - A ) dn 

- M 
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(26) 


?ot  M*/f  Nft  « I , Sq.  (26)  reduces  to 


F = 


M 


*3 


/SyU  <*j  N 


(27) 


i.e()  the  energy  loss  per  cycle  varies  as  the  cute  of  the  maximum  torsional 
norent.  It  is  interesting  to  note  the  olose  correspondence  between  oscillatir. 
torsional  cor.ta.ot  ar.d  oscillating  tangential  contact  of  two  spheres.  In  the 

latter  case,  for  amplitudes  of  tangential  force  small  in  comparison  with 

(?) 

the  energy  loss  per  cycle  is  equal  to 

^ . a - v)  r*3 


36yuafN 


(26) 


is  the  tangential  fowl* 


Fig.  2, 


Fig,  3. 
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